In-plane uniaxial magnetic anisotropy in Ga(Mn)As due to local lattice distortions 

around Mn 2+ ions 



Hemachaiider Subramanian, J. E. Han 

University at Buffalo, State University of New York, Buffalo, NY, 14226 
(Dated: April 13, 2012) 

We theoretically investigate the interplay between local lattice distortions around Mn 2+ impurity 
ion and its magnetization, mediated through spin-orbit coupling of hole. We show that the tetrahe- 
dral symmetry around Mn 2+ ion is spontaneously broken in the paramagnetic regime. Modest local 
lattice distortions around the Mn 2+ ion, along with the growth strain, stabilize magnetization along 
(110) directions, in the ferromagnetic regime. We provide a possible explanation for the experimen- 
tally observed in-plane uniaxial magnetic anisotropy seen in this system by taking into account the 
lack of top-down symmetry in this system. 
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I. INTRODUCTION 

The archetypal dilute magnetic semiconductor (DMS), 
Ga(Mn)As, has been the object of intense investigation 
over the past decade for its technological relevance and is 
also of fundamental theoretical interest ( See [1],[2] and 
references therein). The ferromagnetic and semiconduct- 
ing properties impart multifunctionality to this material 
and it is envisaged to be useful as spintronic material 
in future computational architectures where data stor- 
age and computation could happen at the same place. 
This charge and spin interaction also provides this ma- 
terial with some unusual properties such as temperature 
and hole density dependent magnetic anisotropy 3 - - — . Un- 
der compressive strain imparted by the substrate, this 
material shows both an in-plane uniaxial and cubic mag- 
netic ansiotropy, with the former dominating at high 
temperatures^ - — . This in-plane uniaxial ansiotropy 
has been attributed to surface reconstruction or stack- 
ing faults 1 ^. There is also experimental evidence that 
this uniaxial anisotropy is a bulk property^. In a re- 
cent paper—, Birowska et al suggested that this uniaxial 
anisotropy is due to different nearest neighbor Mn pair 
formation energies along (110) and (110) directions dur- 
ing growth, as observed in their ab initio calculations. Al- 
though there is substantial experimental evidence, there 
are few microscopic theories that provide a convincing 
theoretical explanation for this uniaxial anisotropy. 

In this paper, we investigate the contribution of lo- 
cal lattice distortions around Mn 2+ ion impurity to the 
magnetic anisotropy of Ga(Mn)As. Although, from crys- 
tal field theory perspective, the ground state of Mn 2+ ion 
in Ga(Mn)As is a pure spin state because of its half-filled 
d orbitals, it has been show n 15 i 16 that the d orbitals of 
the Mn 2+ ion hybridize with p orbitals of neighboring 
As ions, and this hybridization should bring the strong 
spin-orbit (SO) effect into the Mn 2+ ion. This opens 
up the possibility of local lattice distortions coupling to 
the spin states of Mn 2+ , resulting in local structural and 
magnetic anisotropy. We expect such coupling with lo- 
cal lattice distortions since it has been observed that the 



spin states of Mn 2+ in Ga(Mn)P and Ga(Mn)N, materi- 
als similar to Ga(Mn)As, couple to local lattice distortion 
modc a 17,18 . Also a photoluminescence experiment 1 ^ on 
Ga(Mn)As reports single-ion uniaxial anisotropy even at 
very low doping concentrations, implying that the origin 
of uniaxial anisotropy lies at the single ion level. 

II. THE MODEL 

In our model, we consider a single Mn 2+ substitutional 
impurity ion bonded to the surrounding four As ions 
placed at the vertices of a tetrahedron, as shown in Fig. 
B^- We assume that the hole electrostatically bound 
to the negatively charged (compared to Ga 3+ sublattice) 
Mn 2+ impurity ion is localized within the M11AS4 subsys- 
tem, and occupies the As p orbitals. The model explicitly 
includes only the triply degenerate Ap orbitals of the four 
As ions and the five degenerate 3d orbitals of the Mn ion. 
These two sets of orbitals are hybridized and this p-d hy- 
bridization is modeled as hopping of electrons between 
As p orbitals and Mn d orbitals^ir— . Hopping parame- 
ters are worked out using the Slater-Koster approach^!. 
Other orbitals are assumed to be the same as in undoped 
GaAs. The Hamiltonian is 

H =H a + H d + H t + £ strain ■ (i) 

The Hamiltonian consists of four parts. Ji a is for holes in 
the acceptor states made of p orbitals of As ions, T~id for 
electrons localized in the d orbitals of the Mn 2+ impurity 
ion, Ht for hopping between the Mn 2+ ion orbitals and 
the acceptor orbitals, and £ strain for local lattice distor- 
tions (treated classically) around the impurity ion. We 
will discuss each of the four terms individually below. 

This model is motivated by similar descriptions of tran- 
sition metal ions in II- VI systems by Ley et ai2&, Mi- 
zokawa and Fujimor i 25 i 29 , and of Mn ion in III-V sys- 
tems by Okabayashi et a l 15 ^ 16 and, Ivanov and Krstajic 
et a P 4 ' 30 . With the exception of Ivanov and Krstajic, 
all others treat the DMS systems as five-ion clusters and 
calculate exchange constants using configuration inter- 
action scheme. These models also assume that the p-d 
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interaction between Mn and the holes in the host lat- 
tice is independent of hole wave-vector. Modifications 
and improvements in the current model over the above 
models are the inclusion of spin-orbit coupling in the ac- 
ceptor part of the Hamiltonian H a , introduction of a lo- 
cal strain energy term to take into account local lattice 
distortions around the Mn 2+ ion, and the use of an intra- 
atomic Coulomb interaction Hamiltonian Hdr^ which is 
rotationally invariant in the spin space. 



being one of A, B, C or D, then the three orbitals are 

\Pi) = \ (\pf) + \pf) + \Pi) + |P?» , i = x,y,z. 

(2) 

These three (six, if spin is taken into account) degener- 
ate p orbitals are assumed to hold one hole in the ground 
state of the H a Hamiltonian. We take the three P or- 
bitals to be the eigenstates of L z operator, which are just 
the Y{~ ,Yi,Yi spherical harmonic functions. 
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FIG. 1: (a) Geometry of the tetrahedral M11AS4. The ions 
denoted as A,B,C and D are As ions. The central Mn 2+ 
has replaced the Ga 3+ ion in Ga(Mn)As. (b) Energy level 
diagram of the MnAs4 subsystem. The electronic configu- 
rations illustrated in the energy level diagram correspond to 
I J = §, Jz = hole spin state and \S = §, S z = §) impu- 
rity spin state. The energy difference between the impurity 
d and the acceptor p states, denoted as E p d, equals e p — e^- 
The dotted lines show all the paths with non-zero hopping 
amplitude (see Table [TTJ) between the impurity ion and the 
acceptor orbitals for the undistorted MnAs4 subsystem. In 
the undistorted case, hopping connects As p orbitals with the 
Mn 2+ d orbitals of t2 g symmetry. 
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In the above equation, cj CT and Ci„ are creation and an- 
nihilation operators for electrons. The second term in- 
troduces spin-orbit coupling to the electrons, where A 
denotes the spin-orbit coupling constant. The subscripts 
j and k in the second term run through the number of 
electrons, and not orbitals. The many-body basis states 
and the eigenstates of H a for the case of five electrons in 
the three P orbitals (corresponding to the presence of a 
single hole) are listed in Appendix A. With the introduc- 
tion of spin-orbit coupling, the six degenerate many-body 
states split into four degenerate J = § states of energy 
5e p — A/2 and doubly degenerate J = \ states of energy 
5e p + A. Thus in our prameterization, the usual spin-orbit 
splitting at the T-point A50 = 3A/2. The value of A, 
derived from A so = 0.341 eV-2 3 ., equals 0.23 eV. We sim- 
ilarly work out the excited many-body eigenstates with 
electron numbers other than five. 



2. Mn 2+ impurity Hamiltonian Hd 



1. Acceptor Hamiltonian rl a 

The first term in the equation fTJ) describes the accep- 
tor states in As ions created due to the substitution of 
a Mn 2+ impurity ion at the position of Ga 3+ ion. The 
hole captured in the hydrogenic acceptor stat o 15 ' 32 is as- 
sumed to reside in orbitals formed by linear combination 
of the atomic p orbitals of the four As ions. At low dop- 
ing concentrations and temperature, it can be reasonably 
assumed that the hole is at k = T-poin^^2^M, 
which implies that the ground state of the hole wavefunc- 
tion has the same phase at all the four As ions. Thus, of 
all possible linear combinations of the twelve atomic As- 
p orbitals (ignoring spin), only three combinations have 
the same phase and will be used in the modeling below. 
If we let \p™) denote the pi orbital of n th As ion, with n 



The impurity ion Hamiltonian describes the intra- 
atomic Coulomb interaction within Mn 2+ d orbitals. 

T-L d =e<2 ^2 n ia + - ^ U Ui^rij-tj 

a 

+ ^2 (U + 2 J) n^mx + -j^2{U - J) n itT n Jcr 

a 

a 

(4) 

In the above equation, d ia , di a and denote creation, 
annihilation and number operators respectively of elec- 
tron with spin a in the z th d orbital of the impurity ion. 
U and J are the intra-ionic direct and exchange Coulomb 
interaction parameters, respectively. The widely used 
Hamiltonian Hd is rotationally invariant in the spin space 
of electrons 3 ^. We choose the ground state occupancy of 
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Parameter 


Value 


Parameter 


Value 


£d — e P 


-4.5 eV 


A 


0.23 eV^ 


U 


1.1 eV-±£ 


J 


0.22 eV±£ 


(pda)o 




(pdn)o 


-0.46 eV-^ 


a 


41.19 N/m2£ 


P 


8.94 N/m^ 



TABLE I: List of parameters used in the model, provided for 
ease of reference. The table lists, by row, energy of d orbitals 
with respect to top of the valence band and spin-orbit cou- 
pling constant, Coulomb direct and exchange integral values, 
Slater-Koster overlap integrals, and, radial and angular strain 
constants used in VFF model, respectively. 



the impurity ion to be five, in accordance with the re- 
sults of numerous experiment s 15 ! 32 ' 34 . The ground eigen- 
state degeneracy of the impurity is six, with the eigen- 
states being the six S z states of the total spin S = |. 
These many-body states are explicitly written down in 
Appendix B. This degeneracy is the result of the spin 
rotational symmetry of Hd- This result is also in accor- 
dance with Hund's first rule of maximizing the total spin 
S of the ion. 

The ground state energy of the Mn 2+ with five elec- 
trons in the d orbitals, denoted by E(d 5 ), is bed + 10U — 
10 J. The center of gravity of the Mn + ion multiplet, 
with six electrons in the d orbitals, denoted by E(d 6 ), is 
6e d + lbU - 5 J, and similarly, E(d 4 ) = 4e d + 6U - 3 J. 
In their photo-excitation experiment^ 5 -, Okabayashi et 
al measured E(d 6 ) + E{d 4 ) - 2E{d b ) (which, in our 
paramctrization, is U + 12 J) to be 3.5 eV (see table 1 of 
the above reference). Using the above value, and assum- 
ing J = 0.2t / 31 ' 35 , we determine the values of U and J 
to be 1.1 eV and 0.22 eV, respectively. Since the ground 
state electronic configuration of Mn ion is foun d 15 ' 32 ' 34 
to be the stable half- filled d 5 configuration, E(d 5 ) must 
be smaller than both E(d 4 ) and E(d 6 ). Thus, with the 
above values of U and J, the value of measured from 
the top of the valence band must be between — 2.7 eV 
and —6.2 eV. We choose — e p to be —4.5 eV, the mean 
value of the above limits. Many-body eigenstates of Hd 
with number of electrons other than five are also worked 
out similarly. Table Q] lists all the parameters used in the 
model, for ease of reference. 



3. Hamiltonian for local lattice distortion, £ s train 

It is theoretically possible for the Mn 2+ impurity to 
induce local lattice distortion s 37 ' 38 . Because of the spin- 
orbit coupling of holes and p-d hybridization between 
holes and the Mn spin, the spin states of the Mn 2+ im- 
purity ion can couple with local lattice distortion a 17 ' 18 ' 39 . 
In this model, we treat the local lattice distortions clas- 
sically We only take into account the distortions of the 
nearest neighbors of Mn 2+ ion, i.e., the four As ions 
A,B,C and D (see Fig. [J). For the five-ion MnAs 4 sub- 
system of Ga(Mn)As, 15 variables are needed to com- 



pletely specify the distorted configuration. Rotational 
and translational invariance remove 6 variables leaving 
9 to completely specify local lattice distortions. These 9 
variables are chosen as follows. Four variables dA, ds, dc 
and do are used to denote the four Mn-As bond lengths, 
leaving five variables to specify angular distortions. We 
choose angles of the bond connecting the Mn and As 
ion B to be fixed, and allow As ion A to move only in 
the A-Mn-B plane; i.e., along the 8 a direction (see Fig. 
[TJa)). The other two ions are allowed to freely rotate 
about the Mn ion in all directions. We have verified that 
choosing a pair of As ions other than A and B for apply- 
ing constraints give the same results (with some choices 
requiring a constraint equation to be solved). 

We parametrize the strain energy of the MnAs 4 
sub-system using Keating's valence force field (VFF) 
model 3 ^: 




E 



3(3d t 



,j = A ) B,C,D 

i>3 



(5) 



Here, d; denotes the position vector of the i As ion, 
and do, the equilibrium Mn-As bond length. The sums 
run over all the four As ions. The parameters a and 
(3 used in VFF model are known for GaAs 3 ^, and are 
41.19 N/m and 8.94 N/m, respectively These values are 
used to evaluate the prefactors in the above equation. We 
take do to be the equilibrium Ga-As bond length. The 
prefactor 2 in front of the radial strain energy term is 
to approximately include the radial strain energy of the 
bonds connecting the As ions A,B,C and D to their near- 
est neighbor Ga ions (next nearest neighbors to Mn 2+ ) 
which arc to be fixed at their equilibrium positions. This 
approximation is equivalent to replacing all the Ga-As 
bonds from a given As ion with a single bond, collincar 
with, and having the same spring constant as that of 
the Mn-As bond. The assumption of a fixed Ga sublat- 
tice around Mn ion is a good approximation, since, an 
X-ray absorption spectroscopy study 4 ^ of the local struc- 
ture around Mn 2+ impurity ion in Ga(Mn)As finds that 
the structural perturbation induced by the impurity ion 
is very localized. With these parameters, the prefactors 
in Eq. © are = 5.7 e V and = 1.3 eV. 

It has to be noted that the distortions couple to the 
electronic degrees of freedom only through the hopping 
amplitudes between d and p orbitals (see Eqs. © and 
© below). In other words, it is assumed that the accep- 
tor Hamiltonian H a is affected by local lattice distortions 
only through its coupling with the Mn 2+ impurity ion. 
This assumption can be rationalized by invoking the dif- 
ference between the strengths of the Mn-As and the Ga- 
As bonds, derived from elastic constants of zinc blende 
MnAs and GaAs. The spring constant a of the Mn- 
As bond, derived from the bulk clastic constants of zinc 
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blende MnAs^i, is about 60 % of the value of spring con- 
stant of Ga-As bond^ 2 -, implying that the Mn-As bonds 
are weaker than the Ga-As bonds. Thus local lattice 
distortions around the Mn + impurity ion can also be 
approximately modeled as the impurity ion itself moving 
with respect to four As ions fixed at their equilibrium 
positions. We verified this by fixing the As ion posi- 
tions and allowing only the Mn 2+ ion to move within 
this fixed four As ion "cage" and reproduced essentially 
the same physics. Thus the assumption that the acceptor 
states are not affected directly by local lattice distortions 
is physically reasonable. Also, with this assumption, we 
achieve considerable computational simplification, by re- 
moving the need for self-consistent calculation of hole 
cigenstatcs in the presence of distortions. This simpler 
strain model, where the Mn 2+ ion moves with respect to 
a fixed As ion cage, provides a simpler way to visualize 
the local lattice distortion modes, as will be shown below. 



4- The hopping Hamiltonian, T-Lt 
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TABLE II: Hopping amplitudes in the absence of local lattice 
distortions. First three d orbitals of t^g symmetry couple with 
atleast one P orbital in the presence of lattice potential Vi a tt, 
whereas the last two d orbitals of e 9 symmetry do not couple 
with any of the P orbitals. 



t + i >yz =—7= X] W^hm k n k (pda) - 2l k m k n k (pdn) 
2 v 2 k=A 

— ^/?>im\n k {pda) — i (l — 2m 2 .) n k (pdir)]. 

(8) 

The hopping amplitude is assumed to depend on bond 
length distortion through modification of hopping param- 
eters (pda) and (pdir) according to Harrison's rule^: 



Finally, the p-d hybridization between As p orbitals 
and Mn 2+ d orbitals is modeled as hopping of electrons 
between the p and d orbitals. 



% — yi u. 



h.c. 
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In Eq. (0, i runs through all the three P orbitals of 
As ion defined in Eq. @, and j runs through the five 
d orbitals of Mn 2+ ion. The hopping amplitudes Uj's 
are functions of the local lattice distortion variables and 
are derived using Slater-Koster formulae 2 ^. To illustrate 
the derivation of these hopping amplitudes, we explicitly 
evaluate i+i,j, z , the hopping amplitude between the Yf 
orbital and yz orbital of Mn 2+ , with the latter assumed 
to be located at the origin of the coordinate system. 



2V2 



^(P+1 


Vlatt 


\dyz) 










c 1 Vlatt 


dyz) " 


i 

-7=2 


(Py\ Vl aU 


\dyz) 






P*\ + 




-{P C X 


+ <tf\) 


Vatt 


\dyz) + 


(7) 


Pf\ + 


(Py\- 


-{P C V 


+ (P?\) 


Vlatt 


\dyz) ■ 





2V2 

If we let (Zfc, m k , n k ) and d k be the direction cosines and 
the length of the bond connecting Mn 2+ ion and k th As 
ion, with k = {A, B, C, D}, we can write the above equa- 
tion for the hopping amplitude t+\^ yz using Slater-Koster 
tables as 



(pda) = (^jj (pda) Q and (pdn) = ^jjj {pdn) . 

(9) 

Here, (pda)o and (pdTr)o arc hopping parameters for 
undistorted Mn-As bonds. We assume (pda)o = 1 eV— , 
and use the relation (pda)o sa — 2.17(p<i7r)(j2i to derive 
(pdTr)o. When the MnAs4 subsystem is undistorted, the 
expression for t + i yz simplfies to 



t 



+i,yz 



V2 



(pda) a _ 2(prf7r) 
3 3V3 



V2 



^- do) 



Similarly 



t- 1>yz =— = and t +1 . yz = -(t-i, yz )*. (11) 
v2 

All the fifteen hopping amplitudes between real orbitals 
for the undistorted MnAs4 subsystem are given in table 
Hi] in terms of to , defined in JTI 



III. 



CALCULATIONAL DETAILS AND 
ANALYSIS 



We solve the Hamiltonian % using second-order de- 
generate perturbation theory, by expanding the ground 
state energy upto second-order in the hopping ampli- 
tudes tij' n 15 ! 25 . The degenerate ground state manifold 
of Wo = Ha + Hd is defined as follows. The ground 
cigenstatcs of the acceptor Hamiltonian H a are the four 
degenerate |j=|,J z = d=|,±|) total angular momen- 
tum eigenstates given in Appendix A. We ignore the 
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J = | split-off band eigenstates. The ground eigen- 
states of the impurity Hamiltonian are the six degen- 
erate \S = §, S z = ±|,±|,±|) eigenstates, given in Ap- 
pendix B. Thus the Hilbert space in which we will work is 
the 24-dimensional product space of the two sets of eigen- 
states mentioned above, i.e., | J = §, J z = ±|, ±i) ® 
|5 = |, 5 Z = ±|, ±|, ±^)- Let S be the energy of these 
24 degenerate states. 

All the terms arising from first-order perturbation ex- 
pansion around Jit vanish trivially, since Ht changes 
the total number of electrons inside the Mn 2+ ion. Wo 
treat paramagnetic case by representing the second- 
order perturbation operator H t ^ Eq L-^ ) % within the 24- 
dimensional Hilbert space spanned by the basis vectors 
mentioned above. We diagonalize the 24 x 24 Hamil- 
tonian matrix to find the ground state eigenvalue, add 
strain energy, and minimize the total ground state en- 
ergy with respect to the distortion variables to find the 
minimum energy configuration of the MnAs4 subsystem. 
We have verified that the Hamiltonian matrix and its 
eigenvectors preserve time-reversal and total spin rota- 
tional symmetry, in the absence of any magnetic field or 
local lattice distortions. 

When the DMS system is ferromagnetic, the Mn 2+ ion 
is magnetized along a particular direction, and we re- 
place its spin operator by a fixed vector. We treat the 
ferromagnetic case by projecting the 24-dimensional ma- 
trix down to a 4-dimensional matrix with only the hole 
spin degrees of freedom included. We then minimize the 
ground state energy ( with strain energy included) with 
respect to distortion variables. Finally, growth- induced 
strain effects are taken into account similarly, by reduc- 
ing the Hilbert space size and working within either light- 
hole or heavy-hole sub-bands, depending upon the type 
of strain. 



Paramagnetic case: 

In the absence of hopping and distortions, the param- 
agnetic ground state degeneracy of the many-body states 
is 24 (6 from 5 = 5/2 states of Mn and 4 from J = 3/2 
states of hole). With hopping, this 24-fold degenerate 
ground state split into 3, 5, 7 and 9-fold degenerate states. 
This splitting can be understood using the total angu- 
lar momentum of the Mn 2+ -hole complex. The quan- 
tum number corresponding to the total angular momen- 
tum F = J + S can take the values 1,2,3 and 4 with 
the degeneracies of 3,5,7 and 9 respectively, and the 
ground state manifold is F = 1, because of the anti- 
ferromagnetic coupling^ between holes and Mn 2+ . Thus 
we observe that the total spin rotational symmetry of 
the entire MnAs4 subsystem is not broken with the in- 
troduction of hopping and the calculated eigenstates are 
also the eigenstates of the total angular momentum op- 
erator F. A heuristic reason for this total spin rotational 
invariance, in the special case of F = 1, is provided in 
Appendix C. This 24-dimensional Hamiltonian matrix, 



whose eigenvectors are also eigenvectors of the total an- 
gular momentum operator F and its z-component F z , is 
then identical, apart from a constant spin-independent 
diagonal term, to the following Hamiltonian involving 
only the spin variables: 

n cfi = Jo J ■ s. (12) 

In the above equation, the components of J (S) refer to 
angular momentum matrices for spins 3/2 (5/2), with 
3/2 (5/2) factored out from the matrices, and Jo is the 
exchange constant. The interaction between Mn 2+ spin 
and the hole spin is anti-fcrromagnctic, and is of Heisen- 
berg type. With our choice of the parameters listed in 
table HI Jo = 0.76 eV— , close to the accepted valued 
of 1.2 ± 0.2 eV. This effective spin Hamiltonian is a di- 
rect consequence of the conservation of the total angular 
momentum F and the second-order perturbation approx- 
imation. 

In the paramagnetic case, the time-reversal symmetry 
of the Hamiltonian is not broken, and we expect the spin 
expectation values of the hole and the Mn 2+ ion to van- 
ish in the ground state, irrespective of distortions. The 
ground eigenstatc energy minimized with respect to the 
local lattice distortions is singly degenerate, and is time- 
reversal symmetric with itself. If we denote the ground 
state eigenvector as |\&gs); then, time-reversal symmetry 
implies 

e|* GS ) = |* GS >. (13) 

Here O is the time-reversal operator. If S denote the spin 
angular momentum operator of cither the hole or Mn + 
ion, then, 

ese^ 1 = -s. (14) 

Therefore, 

(* GS |S|* GS ) =0. (15) 

Thus, even when local lattice distortions are introduced, 
the spin expectation values of the hole and the Mn 2+ ion 
is zero, for the ground eigenstate^l. Distortions break 
orbital, and through spin-orbit coupling, spin rotational 
symmetry, but not time-reversal symmetry. Breaking of 
time-reversal symmetry is needed for the appearance of a 
magnetic dipolc moment. The breaking of spin rotational 
symmetry, due to distortions, is seen at the quadrupolc 
level, with (S" 2 ) ^ (5 2 ) 7^ (S 2 ), in general, for the ground 
cigenstate. 

Energy minimization with respect to distortion vari- 
ables: When we minimize the ground state energy of 
H with respect to all the nine distortion variables, we 
observe spontaneous breaking of tetrahedral symmetry of 
the As ions around Mn 2+ ion. Minimization of energy 
with respect to just radial distortion variables, keeping 
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angles of the Mn-As bonds at equilibrium values, result 
in a configuration where one of the bonds grow shorter 
(5.2% less than the equilibrium value do) than the other 
three bonds (4.8% less than do). The shorter bond can 
be any of the four bonds, and the system has permu- 
tational symmetry with respect to the radial distortion 
variables. The energy landscape has degenerate energy 
minima corresponding to all the four permutations of dis- 
tortions. Figure [5] shows the symmetry breaking when 
the energy is minimized with only two radial strain vari- 
ables, say g?a and de, keeping the other two radial strain 
variables equal and constant, and keeping all angles at 
their equilibrium values. The plot shows that the energy 
is minimal when dA and ds are unequal. 




< d /d / 

FIG. 2: An example of symmetry breaking in MnA&j tetra- 
hedron when the ground state energy is minimized with only 
two radial strain variables, keeping the other two equal and 
constant and all angular variables fixed at their equilibrium 
values. Energy is minimum when dA and da, the length of 
bonds connecting Mn with As ions A and B, are different. Two 
degenerate energy minima, denoted by "X" symbols, seen in 
the plot reflect the permutational symmetry that exists be- 
tween dA and ds. 

Minimization of the ground state energy with only an- 
gular distortion variables result in a broken symmetry 
configuration where three bonds move away from the 
fourth bond (by about 0.4 degrees) while keeping their 
relative angles the same. Again, permutational sym- 
metry holds and the singled-out bond can be any one 
of the four bonds, with the energy surface containing 
four degenerate energy minima corresponding to the four 
possible angular distortion modes. When both radial 
and angular distortion variables are included, the three 
longer bonds move away from the shorter one. Figure [3] 
shows the resultant configuration of the MnAs4 subsys- 
tem. This mode of local lattice distortion is equivalent 
to Mn 2+ ion moving towards the fixed As ion A in the 
simplified strain model introduced above, in which the 
As ions are fixed at their equilibrium positions and the 
Mn 2+ ion is allowed to move within the As ion cage. This 
simplified picture of Mn 2+ ion displacement is approxi- 



mately equivalent to the full five-ion distortion picture. 
With distortion, the triply denegerate F = 1 ground 




FIG. 3: Local lattice distortion when the system is in the 
paramagnetic regime and at one of the four degenerate energy 
minima, exaggerated for clarity. One of the bonds, Mn-A, 
becomes shorter by 0.4% than the other three bonds, while 
the other three bonds move away from the shorter bond by 
0.4 degrees. The tetrahedral symmetry around the Mn 2+ ion 
is spontaneously broken. In the alternate simplified strain 
model in which the only strain variable is the position vector 
of the Mn 2+ impurity ion, the Mn 2+ ion moves towards the 
fixed As ion A. Configurations at other degenerate energy 
minima have the same structure but with the As ion labels 
permuted. 

state manifold splits into a self-time-reversal symmetric 
ground state singlet and an excited state Kramers dou- 
blet. The ground state ceases to be an exact eigenstatc 
of the total angular momentum operator F, and gains a 
small amount of higher total angular momentum, apart 
from F = 1. 

So far, we have not included the growth strain along 
the (001) direction in the model. From experiments, we 
know that in-plane magnetization is the result of com- 
pressive strain exerted by substrates such as GaAs on the 
Ga(Mn)As filmic. This strain splits the J = § mani- 
fold into light-hole (J z = ±4) and heavy-hole (J z = ±|) 
Kramers doublets, with the light-hole doublet remain- 
ing lower in energy 2 -^^. To include the effect of com- 
pressive growth strain, we assume that the holes pop- 
ulate only the light-hole band, and hence include just 
the light-hole J z — ±i doublet in the model, and ne- 
glect the heavy-hole doublet altogether. This results in 
a reduction of the size of our Hilbert space from 24 to 
12. With the inclusion of growth strain, minimization of 
ground state energy with respect to local lattice distor- 
tion result in equal shortening of the four bonds with no 
angular distortions. Thus, symmetry-breaking vanishes 
when only light-holes are considered, in the case of un- 
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polarized Mn ion spin. Realistically, above T c , in the 
bulk DMS sample, the population of holes occupying the 
heavy-hole band will not be negligible, and experiments 
should be able to observe weak local lattice distortions. 
With magnetization of the Mn 2+ ion, distortions reap- 
pear, and are coupled to the direction of magnetization, 
as we will see below. 



Ferromagnetic case: 

At very low temperatures (T <C T c ), the magnetic mo- 
ments of all the Mn 2+ impurity ions align in a particular 
direction, rendering the DMS sample magnetized. We 
model the ferromagnetic case by magnetizing the impu- 
rity ion along a direction parametrized by the polar and 
azimuthal angles (9, $), using rotation operators: 

|X (6, $)) = e -*s.* e -is w e ^ 0j q 5 0j 0) 0] T (16) 

In Eq. (|16p , S y and S z represent angular momentum ma- 
trices for spin | corresponding to rotations about y and 
z-axes respectively, and \x (B, $)) represent the impurity 
ion spinor polarized along (O, <I>) direction. After the re- 
moval of spin degrees of freedom of the Mn 2+ impurity 
ion by magnetizing it, we are left with only the hole de- 
grees of freedom and the size of Hilbcrt space reduces 
further to just four, ignoring the effects of growth strain 
for now. Thus, we diagonalize the 4x4 Hamiltonian ma- 
trix, whose elements are functions of distortion variables 
entering through hopping amplitudes, and impurity ion 
magnetization angles entering through the projection of 
total Hamiltonian onto the S = f spinor \\ (0, $)). We 
then add the radial and angular strain energy E stT ain to 
the ground state eigenvalue to calculate the total ground 
state energy of the system. 

We calculate the structural configuration of the system 
in the ferromagnetic case by minimizing the total ground 
state energy with respect to the local lattice distortion 
variables, for all magnetization angles (0,<I>). Figured] 
shows the total energy of the system as a function of 
the magnetization angles. Eight degenerate minima can 
be observed in the plot, corresponding to the Mn 2+ ion 
magnetization pointing along and opposite to the direc- 
tions of the four As ions. In each of the eight minima, 
the Mn-As bond parallel or anti-parallel to the Mn 2+ 
ion magnetization direction shrinks less (3.3% less than 
do) than the other three bonds (4.7% less than do), and 
the angles between the longer bond and the other three 
bonds is smaller than their equilibrium values, by 1.5 de- 
grees. The distortions in the spin-polarized MnAs4 con- 
figuration move in the opposite direction to that of the 
paramagnetic case, and is shown in Fig. O The min- 
ima corresponding to opposite Mn 2+ ion magnetization 
directions are degenerate and are separated by a finite 
energy barrier. Thus, in the ferromagnetic case with no 
growth strain, magnetic easy axes are along (111) direc- 
tions, when local lattice distortions are included. The 



reason for the appearance of minima when Mn ion 
magnetization points along (111) directions is that, along 
these directions, distortions are more effective in maxi- 
mizing the values of hopping amplitudes between p and d 
orbitals, thereby, decreasing energy. An EXAFS study21 
of Ga(Mn)As has found that the local structure around 
the Mn 2+ ion is significantly altered, and also find that, 
the best fit for their EXAFS oscillations is a model in 
which the Mn 2+ ion is shifted by 0.1 A from the position 
of the Ga ion it has replaced. This experiment corrobo- 
rates our theoretical results. 




$ in units of n 



FIG. 4: Ground state energy of the system as a function of 
magnetization direction of Mn 2+ ion, parametrized by po- 
lar angle <d and azimuthal angle <E>, without the inclusion of 
growth strain. Eight degenerate minima, denoted by "X" 
symbols, corresponding to the magnetization pointing paral- 
lel or anti-parallel to the four Mn-As bonds, i.e., along (111) 
directions, can be seen. The MnAs4 distortion configuration 
corresponding to the left-most minimum at the top is shown 
in Fig. [5] Distortion configurations corresponding to other 
minima have the same structure as in Fig. but with As 
ions A, B, C and D permuted. The lowest energy barrier 
separating any two of these minima is about 2.6 meV. 

To include compressive growth strain, we remove the 
heavy-hole basis states, which are higher in energy than 
light-hole states, from the Hilbert space. The Hilbcrt 
space now has just two dimensions, corresponding to the 
light-hole basis states. The ground state energy, com- 
puted by diagonalizing the 2x2 Hamiltonian matrix, 
with the Mn 2+ ion magnetization pointing along direc- 
tions specified by the angles and $, is minimum when 
the magnetization is in the x-y plane (or = tt/2), i.e. 
perpendicular to the growth direction, in accordance with 
the experiments. We have also verified the case of ten- 
sile strain by including just the heavy-hole doublet, and 
observed uniaxial, out-of-planc magnetization of Mn 2+ 
impurity ion. It has to be noted that the demonstrated 
correspondence between the type of growth strain and 
magnetization direction happens even in the absence of 
local lattice distortions. This correspondence has to do, 
as has been explained also by Sawicki et al&, with the 
maximization of hopping possibilities between Mn 2+ ion 
and the hole orbitals leading to minimization of total 
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FIG. 5: Local lattice distortion when the system is ferromag- 
netic and the magnetization of Mn 2+ ion pointing along Mn-A 
bond direction, without growth strain effects. Distortions are 
exaggerated for clarity. The Mn-A bond becomes longer than 
the other three bonds by about 1.4%, and the angles Mn-A 
bond makes with the other three bonds become smaller by 1.5 
degrees equally. This distortion configuration corresponds to 
the left-most minimum at the top in Fig. 0] The distortion 
in this ferromagnetic case is opposite to that of the paramag- 
netic case shown in Fig. In the alternate simplified strain 
model with fixed As ions, the Mn 2+ moves away from the As 
ion A. 



energy. In the case of tensile strain, the energetically 
low- lying heavy- hole states | J = |, J z = ±|) have more 
hopping possibilities with \S — ^,S Z = T§) than with 
other S z states, thereby lowering the energy of out-of- 
plane polarized spin states of the impurity ion. Sim- 
ilarly, in the case of compressive strain, the light-hole 
states |j=|,J z = ±^) have more hopping possibilities 
with \S = \,S Z = T5) than with other S z states, and 
hence lower the energy of in-plane spin-polarized states 
of the impurity ion. Restricting the Hilbert space to just 
the light-hole doublet results in magnetization pointing 
in the x-y plane, all directions in the plane being equally 
possible. Local lattice distortions further reduce the in- 
plane rotational symmetry and lead to bi-axial magnetic 
anisotropy. Figure [5] shows the total energy of the system 
as a function of magnetization angles and <f>. As can be 
seen in the figure, there are four minima corresponding 
to the magnetization of Mn 2+ ion pointing along (110) 
directions. Variation of energy when the magnetization 
is rotated in the x-y plane is shown in Fig. [71 Local lat- 
tice distortion at one of these minima is shown in Fig. [8] 
The figure shows local lattice distortions when magneti- 
zation of the Mn 2+ ion points along (110). In this case, 
the bonds connecting As ions C and D are shorter (4.1% 
less than do) than that of A and B (3.5% less than do), 
and the angle between C and D increases above cos _1 (^), 




<i> in units of n 



FIG. 6: Ground state energy of the MnAs4 subsystem as 
a function of the direction of magnetization of the Mn 2+ 
ion, parametrized by O and when the system is ferro- 
magnetic and has growth strain that splits heavy and light 
holes. With the inclusion of only the light holes, minimum 
ground state energy occurs along (110) directions (6 = 7r/2, 
$ = 7r/4, 37r/4, 57r/4, 77r/4), as can be seen in the figure, and 
are denoted by "X" symbols. Local lattice distortion configu- 
ration corresponding to the first minimum from left is shown 
in Fig. [S] Distortion configurations corresponding to other 
minima have the same structure as in Fig. but with As 
ions A, B, C and D permuted. 




$ in units of k 



FIG. 7: Energy as a function of magnetization angle $ of 
the Mn 2+ ion in the x-y plane, in a sample with compres- 
sive growth strain. The minima correspond to (110) mag- 
netization directions (0 = n/2, $ = 7r/4, 37r/4, 57r/4, 77t/4). 
The system shows biaxial anisotropy, with the introduction 
of strain-induced heavy hole - light hole splitting. The energy 
barrier between two minima is 4.5 meV. 

whereas, the angle between A and B decreases, both by 
about 1.6 degrees. In the case of (110) magnetization, 
the roles of A and B are replaced by C and D, and vice 
versa. The Mn 2+ ion moves towards the x-y plane in the 
equivalent simplified strain model where all other ions 
are fixed at their equilibrium positions. 

We reproduce the experimentally observed uniaxial 
anisotropy if we include the gradual change in the bond 
lengths of Ga(Mn)As along (001) growth direction, from 



/ 

FIG. 8: Local lattice distortion when magnetization of the 
Mn 2+ ion points along (110) direction. All bonds decrease in 
length, with bonds connecting C and D ions decreasing more 
(4.1% less than do). The angles between C and D increases, 
whereas, the angle between A and B decreases, both by 1.6 
degrees. Distortions are exaggerated for clarity. The Mn 2+ 
ion moves towards the x-y plane in the alternate strain model 
with fixed As ions. 



the substrate to the surface. To take this gradual change 
of bond lengths into account, in the case of compressive 
strain, we set the equilibrium bond lengths of As ions C 
and D shorter than that of As ions A and B by 1%, as 
in Fig. [SI and minimize the energy with respect to dis- 
tortion variables. Figure [5] shows the variation in total 
ground state energy as the magnetization of Mn 2+ ion 
is rotated in the x-y plane. We observe the breaking of 
the remaining bi-axial symmetry, with the system favor- 
ing (110) direction of magnetization over (110) direction. 
Zemen et a& refer to an experiment where Ga(Mn)As 
grown on top and bottom surfaces of the same GaAs sub- 
strate show perpendicular in-plane magnetizations. This 
observation can be simply explained within this model 
once we realize that, for Ga(Mn)As grown on the top sur- 
face of GaAs, As ion bonds C and D are shorter, whereas, 
for Ga(Mn)As grown on the bottom surface, bonds A and 
B are shorter. 



IV. SUMMARY AND CONCLUSIONS 

In this paper, we have built a many-body tight-binding 
model for Ga(Mn)As DMS system to explain the exper- 
imentally observed uniaxial anisotropy. The model in- 
cluded spin-orbit interaction, intra-atomic Coulomb in- 
teraction inside the Mn 2+ impurity ion, hopping between 
the impurity ion and hole states and local lattice distor- 
tions around the impurity ion. The theory applies in the 
low doping limit (near k = 0, T-point) where we assumed 
symmetric wavefunction of As p orbitals. We have com- 



FIG. 9: Ground state energy as a function of magnetization 
angle of the Mn 2+ ion in the x-y plane, when compressive 
growth strain is included, both at the level of band struc- 
ture by choosing light-holes, and at the level of local lattice 
structure around the impurity by choosing different equilib- 
rium bond lengths for As ions in the +z and — z directions. 
The blue curve shows energy variation when the equilibrium 
lengths of bonds A and B are chosen to be shorter than that 
of C and D. The system shows uniaxial magnetic anisotropy 
along (110) direction. The red curve shows the same scenario 
with C and D shorter than A and B, resulting in anisotropy 
along (110) direction. 



putationally derived the Heisenberg exchange interaction 
between the Mn 2+ impurity ion and the hole spins in the 
distortionless case. Spontaneous tetrahedral symmetry 
breaking around the impurity ion is observed, even when 
the system is paramagnetic, when strain energy due to 
local lattice distortions are included in the Hamiltonian. 
Magnetization of impurity ion couples to the local lattice 
distortions through spin-orbit coupling to produce the ex- 
perimentally observed anisotropies. Under compressive 
(tensile) strain, we observe in-plane (out-of-plane) mag- 
netic anisotropy, by confining ourselves to light (heavy) 
holes. Growth strain, or, in other words, the breaking of 
top-down symmetry, also leads to the breaking of in-plane 
biaxial symmetry and results in uniaxial anisotropy. 

There is considerable scope for this model to be ex- 
panded further, through inclusion of non-zero k states of 
holes, non-zero temperatures and other ionization states 
of Mn 2+ ior>££, to capture the full complexity of the very 
interesting anisotropy behavior of Ga(Mn)As, as a func- 
tion of temperature, hole density and Mn concentration, 
apart from the strain from substrate. 
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V. APPENDIX 

Basis and eigenstates of the acceptor 
Hamiltonian H a 



Negative S z states are obtained by just nipping the spins 
in the many-body states from up to down and vice versa. 



In the absence of hopping, spin-orbit coupling and fer- 
romagnetism, the six degenerate many-body basis states 
of H a are denoted as 



in, n, t), in, u, i>, in, i, u) 
in, t, u), n, u,u), it, u,u) 



(17) 



with the three P orbitals being just the Y{~ ,Yi,Yi~ 
spherical harmonic functions. We work out the eigen- 
states of the hole in the presence of spin-orbit coupling, 
by diagonalizing H a in the basis of Eq. (fT7| . 



J=- 2 ,J, 

J=-,J, 
2 

J=-,J» 
2 

2 



J=- 2 ,J Z 



-3 
-1 

i 

2 
3 
2 

T 
l 

2 



= ltUtU>, 

= ^(m,n,t>-V2|tj,4.n>), 

-^(i4,n,n)-V2in,t,n) 

it, n,ti), 
i 



v/3 



v%4.,t4.,t> + ltU,t4.) 



l 

71 



^/2|;,tl,U) + IU,t,U> 



(18) 



B. Eigenstates of the Mn 2+ impurity Hamiltonian 

H d 

The ground eigenstates of Hd, when five electrons oc- 
cupy the d orbitals, are given below. The five d orbitals 
in the following many-body states are xy, yz, zx, x 2 — y 2 
and z 2 . 



2' 2 



lt,t,t,t,t> 



(19) 



S 



,s x = 



s 



:,S Z 



V5 



[lt,t,t,t,l) + lt,t,t,l,t) + 



lt,t,i,t,t) + lt,l,t,t,t) + ll,t,t,t,t)] 

(20) 

[|t,t,t,t,t) + lt,t,l,t,l) + 



|t,t,t,t,t) + lt,t,t,t,t) + lt,t,t,t,t) 
|t,t,t,t,t) + lt,t,t,t,t) + lt,t,t,t,t) 
||,t,|,t,t) + ll,t,t,t,t)] 

(21) 



A heuristic understanding of the effective spin 
Hamiltonian Eq. (Tl2|) : 



To motivate the total spin rotational invariant form of 
the effective Hamiltonian Eq. (fT2|) . we write down the 
three eigenstates of the total angular momentum opera- 
tor F, \F = 1,F Z = 0,±1), in terms of \S Z ,J Z ): 

if=i,F,=+i>=-vT|i.^)+vTf|i.T>-V9&lfi)+v^lT.«>. 
|F=i A =o>=yi|| ) ^)-^x|i,^)+y*|^,j)-^|:|^,|), 

\ F = 1 ' F ^=^ 1 )=Vw \ 5'Tr) _ \/5'l T i ^T i )+V^rol"2^4)+\/¥| tt'I, 

(22) 

To prove that these three eigenstates of F are also the 
three lowest eigenstates of the full Hamiltonian H in the 
undistorted case, we first note that the off-diagonal ele- 
ments oiHt (Eg-no) ^ 1 between the three eigenstates (f2"2")l 
are zero. This follows from the structure of the hopping 
amplitudes given in table[Hl which shows that the double- 
hopping process conserves the z-component of the total 
angular momentum, F z = S z + J z . Double- hopping pro- 
cess, irrespective of local lattice distortions, always con- 
serve the z-component of the total spin angular momen- 
tum, i.e., AS Z + As z ole = 0. In the undistorted case, in 
addition, z-component of the orbital angular momentum 
of the final hole state is the same, after double-hopping, 
as that of the initial hole state, because the P orbital 
involved in both the hop-out and hop-in processes is the 
same (please see table|TT]), and so, AL\ ole = 0. Therefore, 
AF Z = AS Z + As h z ole + AL h ° le = AS Z + AJ Z = 0. This 
proves that the off-diagonal elements of the second-order 
perturbation operator between different F z eigenstates 
are zero. 



To prove that the three diagonal matrix elements of 
Ht rg~^j % between F z — 0, ±1 are degenerate, we in- 
voke the rotational symmetry of the ground eigenstates 
of H a and Hd, and the fact that the fifteen hopping am- 
plitudes appear in the Hamiltonian matrix elements only 
as 



— xy,yz,zx, 

2 2 2 

x —y ,z 



(23) 



In the undistorted case, only the three double hopping 
amplitudes, Too, T +1 +1 and T_i _x, are non-zero, and 
are all equal to t^. This results in the three diagonal 
matrix elements to be the same. Thus, the operator 
Ht (Ed-Ha) ^ 1 IS a constant times identity matrix in the 
sub-space of total angular momentum F = 1 , in the case 
of undistorted MnAs4 subsystem. 



11 



1 Igor Zutic, Jaroslav Fabian, S. Das Sarma, Rev. Mod. 
Phys. 76, 323 (2004) 

2 T. Jungwirth, Jairo Sinova, J. Masek, J. Kucera, and A. 
H. MacDonald, Rev. Mod. Phys. 78, 809 (2006) 

3 J. Zemen, J. Kucera, K. Olejnik, and T. Jungwirth, Phys. 
Rev. B 80, 155203 (2009) 

4 Priya Mahadevan, Alex Zunger, and D. D. Sarma, Phys. 
Rev. Lett. 93, 177201 (2004) 

5 K.Y. Wang, M. Sawicki, K. W. Edmonds, R. P. Campion, 
S. Maat, C. T. Foxon, B. L. Gallagher and T. Dietl, Phys. 
Rev. Lett. 95, 217204 (2005) 

6 C. Timm and A. H. MacDonald, Phys. Rev. B 71, 155206 
(2005) 

7 T. Dietl, H. Ohno and F. Matsukura, Phy. Rev. B 63, 
195205 (2001) 

8 M. Sawicki, F. Matsukura, A. Idziaszek, T. Dietl, G. M. 
Schott, C. Ruester, C. Gould, G. Karczewski, G. Schmidt, 
and L. W. Molenkamp, Phys. Rev. B 70, 245325 (2004) 

9 M. Sawicki, K.Y. Wang, K. W. Edmonds, R. P. Campion, 
C. R. Staddon, N. R. S. Farley, C. T. Foxon, E. Papis, E. 
Kamiaska, A. Piotrowska, T. Dietl and B. L. Gallagher, 
Phy. Rev. B 71, 121302 (2005) 

10 U. Welp, V. K. Vlasko-Vlasov, X. Liu, J. K. Furdyna and 
T. Wojtowicz, Phys. Rev. Lett. 90, 167206 (2003) 

11 S. C. Masmanidis, H. X. Tang, E. B. Myers, Mo Li, K. 
De Greve, G. Vermeulen, W. Van Roy, and M. L. Roukes, 
Phys. Rev. Lett. 95, 187206 (2005) 

12 Unjong Yu, Abdol-Madjid Nili, Karlis Mikelsons, Brian 
Moritz, Juana Moreno and Mark Jarrell, Phys. Rev. Lett. 
104, 037201 (2010) 

13 M. Kopecky, J. Kub, F. Maca, J. Masek, O. Pacherova, 
B. L. Gallagher, R. P. Campion, V. Novak, T. Jungwirth, 
Phys. Rev. B 83, 235324 (2011) 

14 M. Birowska, C .Sliwa, J. A. Majewski and T. Dietl, 
larXiv: 1202.32951 ^1 [cond-mat.mtrl-sci] 

15 J. Okabayashi, A. Kimura, O. Rader, T. Mizokawa, A. 
Fujimori, T. Hayashi and M. Tanaka, Phys. Rev. B 58, 
4211 (1998) 

16 J. Okabayashi, A. Kimura, T. Mizokawa, A. Fujimori, T. 
Hayashi and M. Tanaka, Phys. Rev. B 59, 2486 (1999) 

17 G. Hofmann, F. G. Anderson, J. Weber, Phys. Rev. B 43, 
9711 (1991) 

18 A. Stroppa and G. Kresse, Phys. Rev. B 79, 201201, (2009) 

19 O. M. Fedorych, E. M. Hankiewicz, and Z. Wilamowski, 
Phys. Rev. B 66, 045201 (2002) 

20 Figures created using the program "Balls & Sticks". De- 
tails can be found in the paper - T. C. Ozawa, Sung J. 
Kang, Balls & Sticks: Easy-to-Use Structure Visualization 
and Animation Creating Program J. Appl. Cryst. 37, 679 
(2004) 

21 F. D. M. Haldane and P. W. Anderson, Phys. Rev. B 13, 
2553 (1976) 

22 B. E. Larson, K. C. Haas, H. Ehrenreich and A. E. Carls- 
son, Phys. Rev. B 37, 4137 (1988) 

23 A. K. Bhattacharjee, Phys. Rev. B 46, 5266 (1992) 

24 P. M. Krstajic, F. M. Peeters, V. A. Ivanov, V. Fleurov, 
K. Kikoin, Phys. Rev. B 70, 195215 (2004) 

25 T. Mizokawa and A. Fujimori, Phys. Rev. B 56, 6669, 
(1997) 



2b T. O. Strandberg, C. M. Canali, A. H. MacDonald Phys. 
Rev. B 81, 054401 (2010) 

27 J. C. Slater, G. F. Koster, Phys. Rev. 94, 1498 (1954) 

28 L. Ley, M. Taniguchi, J. Cxhijsen, R. L. Johnson and A. 
Fujimori, Phys. Rev. B 35 2839, (1987) 

29 T. Mizokawa and A. Fujimori, Phys. Rev. B 48, 14150, 
(1993) 

30 V.A. Ivanov, P.M. Krstajic, F.M. Peeters, V. Fleurov and 
K. Kikoin, J. Mag. Mag. Mater. 258, 237, (2003) 

31 A. M. Oles, Phys. Rev. B 28, 327 (1983) 

32 M. Linnarsson, E. Janzten, B. Monemar, M. Kleverman, 
and A. Thilderkvist, Phys. Rev. B 55, 6938 (1997) 

33 M. Cardona, N. E. Christensen and G. Fasol, Phys. Rev. 
B 38, 1808, (1988) 

34 J. Szczytko, A. Twardowski, K. Swiatek, M. Palczewska, 
M. Tanaka, T. Hayashi and K. Ando, Phys. Rev. B 60, 
8304 (1999) 

35 J. Friedel and C. M. Sayers, J. Phys. (Paris) 38, 697, (1977) 

36 Jose Luis Martins and Alex Zunger, Phys. Rev. B 30, 6217 
(1984) 

37 I. N. Demchenko, K. Lawniczak-Jablonska, T. Story, V. 
Osinniy, R. Jakiela, J. Z. Domagala, J. Sadowski, M. 
Klepka, A. Wolska and M. Chernyshova, J. Phys.: Con- 
dens. Matter 19, 496205 (2007) 

38 S. Katsumoto et al, Mater. Sci. Eng. B, 84, 88 (2001) 

39 X. G. Guo, X. S. Chen, Y. L. Sun, X. H. Zhou, L. Z. Sun, 
J. C. Cao and W. Lu, Phys. Rev. B 69, 085206 (2004) 

40 P. N. Keating, Phys. Rev. 145, 637 (1966) 

41 F. d'Acapito, G. Smolentsev, F. Boscherini, M. Piccin, G. 
Bais, S. Rubini, F. Martelli, and A. Franciosi, Phys. Rev. 
B 73, 035314 (2006) 

42 J. Qi, J. A. Yan, H. Park, A. Steigerwald, Y. Xu, S. N. 
Gilbert, X. Liu, J. K. Furdyna, S. T. Pantelides and N. 
Tolk, Phys. Rev. B 81, 115208 (2010) 

43 Walter A. Harrison, Elementary Electronic Structure, 
World Scientific, 1999 

44 Walter A. Harrison, Electronic Structure and the Proper- 
ties of Solids, Dover Publications, 1988, p-451 

45 H. Ohno, J. Magn. Magn. Mater. 200, 110 (1999) 

46 The reason that the value of Jo is different from the ac- 
cepted value is because of our inclusion of spin-orbit cou- 
pling term in the Hamiltonian, whereas it is absent in the 
model of okabayashi et al^, from which we derive the val- 
ues for our parameters. 

47 The fact that the expectation value of Mn 2+ spin van- 
ishes under time-reversal symmetric condition makes the 
assumption of Mn 2+ spin as a classical vector with a fi- 
nite magnetic moment untenable under the same con- 
dition. Thus this assumption cannot be made in non- 
ferromagnetic cases or when a system is frustrated, with 
zero net magnetization. 

48 A. Shen, H. Ohno, F. Matsukura, Y. Sugawara, N. Akiba, 
T. Kuroiwa, A. Oiwa, A. Endo, S. Katsumoto, and Y. Iye, 
J. Cryst. Growth 175, 1069 (1997) 

49 X. Liu, Y. Sasaki, and J. K. Furdyna, Phys. Rev. B 67, 
205204 (2003) 

50 V.F. Sapega, M. Moreno, M. Ramsteiner, L. Dalweritz, 
and K. Ploog, Phys. Rev. B 66, 075217 (2002) 



